Double resonance frequency shift in a hydrogen maser 



M.A. Humphrey, D.F. PhiUips and R.L. Walsworth 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 
(February 2, 2008) 



We use the dressed atom formalism to calculate the frequency shift in a hydrogen maser induced 

' by applied radiation near the Zeeman frequency, and find excellent agreement with a previous 

\ calculation made in the bare atom basis. The maser oscillates on the = f , Amp = hyperfine 

■ transition, while the applied field is swept through the F = 1, Am^ = ±1 Zeeman resonance. We 

' determine the effect of the applied field on the Zeeman levels using the dressed atom picture, and 

1— ( , then calculate the maser frequency shift by coupling the dressed states to the microwave cavity. 

, ^ ■ Qualitatively, the dressed-atom analysis gives a new and simpler physical interpretation of this 

' ' ' double resonance process, which has applications in precision hydrogen Zeeman spectroscopy, e.g., 

. in fundamental symmetry tests. 
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I. INTRODUCTION 



. Since its development nearly 40 years ago, the hydrogen maser has served as a robust tool capable of making high 
^ " precision measurements by utilizing its excellent frequency stability [|l|-^. Hydrogen masers are currently the most 
^ • stable active oscillators over intervals of seconds to days, with applications including very long baseline interferometry, 
' ^ ' deep space tracking and navigation, and metrology. In addition, hydrogen masers have been used to make precision 
J>-j. atomic physics measurements [^|-|8| and for sensitive tests of general relativity |^ and quantum mechanics [ p^p| ] . 

■ A hydrogen maser operates on the first-order magnetic field- independent AF — 1, Arn^? = hyperfine transition, 
between states |2) and |4) of the electronic ground state (see Fig. Q). As shown in Fig. ^ atoms in states |f) and |2) are 
injected into a storage bulb residing in a microwave cavity tuned near the hyperfine frequency, creating the population 
inversion necessary for active maser oscillation. The microwave field stimulates a small, coherent magnetization in 
the atomic ensemble, and this magnetization acts as a source to increase the microwave field. With sufficiently high 
atomic flux and low cavity losses, this feedback induces active maser oscillation. The maser frequency is very stable; 

' a well engineered hydrogen maser can have fractional frequency stability on the order of lxfO~^^ over intervals of a 
: few hours. 

' Hydrogen masers can also be used as sensitive probes of the F — 1, Amp = ±f Zeeman resonances through 

■ a double resonance technique jl^ , in which an oscillating magnetic field tuned near the atomic Zeeman resonance 
shifts the AF = 1, Amp = maser frequency. At low static magnetic fields, this maser frequency shift is an 



^ ■ antisymmetric function of the detuning of the applied oscillating field from the Zeeman resonance. Thus, by observing 
the antisymmetric pulling of the otherwise stable maser frequency, the hydrogen Zeeman frequency can be determined 



with high precision. 

^ '' An early investigation of atomic double resonance was made by Ramsey ]l3[ , who calculated the frequency shift 
O |. between two levels coupled by radiation to other levels. This calculation treated the problem perturbatively to first 
■ order in the coupling field strength, and it neglected damping. Later, Andresen |l^,|f^ calculated the frequency shift 
• I— I , in a hydrogen maser due to an applied field oscillating near the F = f Zeeman frequency (see Sec. and then 
' measured this double resonance effect using a hydrogen maser, finding reasonable agreement with the calculation. 
Andresen's calculation employed a bare atom basis, treating the problem to second order in the applied field strength 
and including phenomenological damping terms, but used an oversimplified description of spin-exchange relaxation. 
Savard [|l5| revisited the problem with a more realistic spin-exchange relaxation description jl^, and found a small 
correction to the earlier work. 

Although the work of Andresen and Savard provides a complete description for the double resonance maser frequency 
shift, intuitive understanding is obscured by the length of the calculations and the use of the bare atom basis. 
In particular, these works demonstrate that the amplitude of the antisymmetric maser frequency shift is directly 
proportional to the electronic polarization of the masing atomic ensemble. The maser frequency shift vanishes as this 
polarization goes to zero. The previous bare atom analyses provide no physical interpretation of ths effect. 
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Since the dressed atom formalism [|7| often adds physical insight to the understanding of the interaction of matter 
and radiation, we apply it here to the double resonance frequency shift in a hydrogen maser. In a two step process, we 
first use the dressed atom picture to determine the effect of the applied Zeeman radiation on the atomic states. Then, 
we analyze the effect of the microwave cavity field on the dressed states and determine the maser frequency shift. 
We find excellent quantitative agreement between the dressed atom calculation and the previous bare atom result. 
We exploit the dressed atom calculation to offer a simple physical interpretation of the double resonance hydrogen 
maser frequency shift, including a straightforward explanation of the polarization dependence. We conclude by noting 
the experimental application of the double resonance effect to perform improved, high precision hydrogen Zccman 
spectroscopy. 



II. MASER OSCILLATION FREQUENCY 

In a hydrogen maser (Fig. molecular hydrogen is dissociated in an rf discharge and a thermal beam of hydrogen 
atoms is formed. A hexapole state selecting magnet focuses the higher energy, low-field-seeking hyperfine states (|1) 
and |2)) into a quartz maser bulb at about 10^^ atoms/sec. Inside the bulb (volume « 10^ cm^), the atoms travel 
ballistically for about 1 second before escaping, making « 10'' collisions with the bulb wall. A Teflon coating reduces 
the atom-wall interaction and thus inhibits decoherence of the masing atomic ensemble by wall collisions. The maser 
bulb is centered inside a cylindrical TEqh microwave cavity resonant with the 1420 MHz hyperfine transition. As 
described above, the microwave field and atomic ensemble form a coupled system with positive feedback that will 
actively oscillate near the AF = 1, Amp — hyperfine frequency if there is a sufficiently high atomic flux into the 
maser bulb. Since the atoms are confined to a region of uniform microwave field phase, their velocity is effectively 
averaged to zero over the interaction time with the microwave field, and first-order Doppler effects are eliminated. The 
maser signal is inductively coupled out of the microwave cavity and amplified with an external receiver. Surrounding 
the cavity, a solenoid produces a weak static magnetic field to maintain the quantization axis inside the maser bulb, and 
a pair of Helmholtz coils produces the oscillating transverse magnetic field that drives the _F = 1 Zeeman transitions. 
The cavity, solenoid, and Zeeman drive coils are all enclosed within several layers of magnetic shielding. 

To calculate the maser oscillation frequency, we first determine the relationship between the microwave field and 
a source magnetization. If we assume that a single cavity mode is dominant within the microwave cavity, we may 
write the microwave magnetic field as H(r,cj) — \/47rpc(w)Hc(r), where pc characterizes the frequency- (i.e. time- 
) dependent amplitude and He represents the time-independent spatial variation of the mode. A straightforward 
application of Maxwell's equations relates pc{^) to a source magnetization M(a;): 

{^l + ^ - ^')Pc(^) - V^Lo\H\^^),VkM{uj) (1) 
Qc 

(z) 

where ui is the maser frequency, luq is the cavity frequency, {H(-. )t is the average of the microwave field's z-component 
over the bulb, and Vb is the volume of the maser bulb. The second term on the left has been introduced phenomeno- 
logically to account for cavity losses characterized by Qc, the cavity quality factor. The source magnetization is 
produced by the atomic ensemble, and is given by the expectation value of the magnetic dipole operator (neglecting 
the term precessing at —uj), 

M{lu) =N {fi)^N Tr(pii) = N^i2APA2{oj) (2) 

where N is the atomic density and /i24 is the hyperfine transition dipole matrix clement, /^24 ~ ^fJ-B- The atomic 
coherence pi2(j-^) is found by solving the Bloch equations 

P= ^[P, Hq] + - [p, Hint] + -f Prelax (3) 

in steady state. The unperturbed Hamiltonian, Hq, includes only the atomic state energies (given by the hyperfine 
and Zeeman interactions), and the interaction Hamiltonian, Hint, describes the effect of the microwave cavity field, 
which couples states |2) and |4). For perfect hexapole state selection, the flux term is written 

P/;«. = ^(|l)(l| + |2)(2|)-rp, (4) 

where the positive terms account for the injection of atoms in states |1) and |2) at rate r, and the last term accounts 
for bulb escape. The relaxation term, Preiax, phenomenologically describes population and coherence relaxation due 
to wall collisions, magnetic field inhomogeneities, and spin-exchange collisions. 
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In the absence of Zeeman radiation and for small detuning of the cavity from the |2) ^ |4) hyperfine transition 
frequency, 1^24, the maser frequency is found using Eqns. - 1 to be 

Lj = CJ24 + ^ i^^c - ^24) (5) 

where Qi — 1JJ24/AW, the atomic line-Q, is on the order of 10^ Since the cavity-Q is on the order of 10'^, the cavity 
pulling is highly suppressed, and the maser oscillates near the atomic hyperfine frequency, a;24- 

In the presence of Zeeman radiation, the maser frequency is shifted, as first shown by Andresen Jl2| and measured 
anew here (see Fig. His calculation of the double resonance shift included the applied Zeeman field in the 

interaction Hamiltonian but otherwise left unchanged the above analysis for the maser frequency. To second order in 
the Rabi frequency of the applied Zeeman field, IX12I, and in terms of the unperturbed maser Rabi frequency l-'^^jl, 
atom flow rate r, population decay rate 71, hyperfine decoherence rate 72, and Zeeman decoherence rate 7^, the small 
static field limit of the maser shift is given by 

where 5 is the detuning of the applied field from the atomic Zeeman frequency, and pJi ~ P33 — ^/ (271) is the steady 
state population difference between states |1) and |3) in the absence of an applied Zeeman field (following Eqn. 8 of 
ref. Physically, the population difference between states |1) and |3) represents the electronic polarization of the 

hydrogen ensemble p^]: 

P=^-^=2Tr{pSz) = Pii-pz^. (7) 

Equation]^ implies that a steady state electronic polarization, and hence a population difference between states |1) 
and 1 3) injected into the maser bulb, is a necessary condition for the maser to exhibit a double resonance frequency 
shift. Walsworth et. al. demonstrated this polarization dependence experimentally by operating a hydrogen maser in 
three configurations: (i) with the usual input flux of atoms in states |1) and |2); (ii) with a pure input flux of atoms 
in state |2), where the maser frequency shift vanishes; and (iii) with an input beam of atoms in states |2) and |3), 
where the maser shift is inverted [Q. 

For typical applied Zeeman field strengths, the 1420 MHz maser frequency is shifted tens of mHz (sec Fig. ||), a 
fractional shift of « 10^^^. However, the shift is easily resolved because of the excellent fractional maser frequency 
stability (parts in lO^'^). 



III. DRESSED ATOM CALCULATION 



We now consider an alternative approach to calculating the double resonance maser frequency shift using the dressed 
atom picture. We retain the Maxwell-Bloch formalism of Eqns. - H; however we determine the steady state coherence 
^42(0;) in a dressed atom basis, including the atomic state energies, the applied Zeeman field, and the atom/Zeeman 
field interaction. For simplicity, we assume the static magnetic field is sufficiently low that the two -F = 1, Amj? = ±1 
Zeeman frequencies are nearly degenerate, W12 — W23 ^ 7z, as is the case for typical hydrogen maser operation. We 
use the simplified spin-exchange relaxation model and neglect the small spin-exchange correction to the double 
resonance maser frequency shift [p^ . 



A. Dressed atom basis 



By incorporating the applied Zeeman field into the unperturbed Hamiltonian, it takes the form Hq = Ha + Hf + Vaf. 
The atomic states (defining state |2) as energy zero) are described by Ha = ?iaji2|l)(l| — ?iw23|3)(3| — ?iu;24|4)(4|; the 
applied Zeeman field (at frequency ujz) is described by Hf = TiuzO'^O'', and the interaction between them is given by 

Vaf=ng{a + a^) (2| + |2) (3| + /i.e.] . (8) 

Here, the Zeeman field creation and annihilation operators are and a, h.c. denotes Hermitian conjugate, and g is 
the single-photon Rabi frequency for the Zeeman transitions. We will use eigenkets with two indices to account for 
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the atomic state and the number of photons in the Zeeman field, denoted by n. We select four of these as our bare 
atom/Zeeman field basis, {|1, n — 1), |2, n), |3, n + 1), |4, n)}, where the first entry indicates the atomic state and the 
second entry indicates the Zeeman photon number. We note that for a resonant field, cuz = i^i2, the first three basis 
states are degenerate. Also, n 3> 1 in practice for there to be any measureable double resonance maser frequency 
shift. 

In this bare atom/Zeeman field basis, the unperturbed Hamiltonian operator takes the following matrix form: 
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(9) 



where 8 — ujz — i^i2 is the detuning of the applied Zeeman field, and 

X\2 — / ^23 

— — g\/n « g^/n + 1 = — — 



(10) 



define the Zeeman field Rabi frequency (the factor of two has been inserted to be consistent with our rotating wave 
approximation convention; see note p8[). By diagonalizing iJo, we find new basis states which physically represent 
the atomic states dressed by the applied Zeeman field. The dressed energy levels are the eigenvalues of Hq: Ea = h^, 

Ell — 0, Ec = —hft, and E4 — —hu}2A, where 51 = \J 5'^ + 5^12 represents the generalized Rabi frequency. The dressed 
states are the eigenvectors of Hq: 
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\b) 
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(11) 



1 3), while in the limit of large positive 5, 



Note that in the limit of large negative (5, |a) |1) and |c 
|c) — > This will become important in the physical interpretation of the maser frequency shift (see Sec. fV). 
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An operator O transforms between bare atom to dressed atom bases as = T~^0''T, where T is the unitary 
matrix linking the dressed atom and bare atom basis states (coefficients of Eqn. |l^). The dressed and bare atom 
energies and eigenvectors are equivalent for the F = Q hyperfine state |4) because this state is unaffected by the 
applied Zeeman field. 



B. Dressed basis Bloch equations 

We now couple the dressed states to the microwave cavity using the Bloch equations, which remain of the form 

={[P\ H^] + J Htt] + P\ela. + P^flu. ■ (12) 

The unperturbed Hamiltonian, iJg, now accounts for the bare atom energies and the applied Zeeman driving field, 
while the microwave cavity field is included in the interaction Hamiltonian, Hint- Since the dressed states |a), \b), and 
|c) all have a component of the atomic state |2) (see Eqn. ^l]), the microwave field couples state |4) to each: 

Htt = §f ^24|a)(4| + ^H2,\b){A\ - §f i^24|c)(4| + h.c. (13) 

Note that H24, = (2|/i • Hc|4) is the only nonzero coupling between bare atom states supported by the TEqh mode 
microwave cavity. To transform the relaxation terms into the dressed basis, we make the approximation that all 
relaxation rates (population decay 71, hyperfine decoherence 72, and Zeeman decoherence 7^) have the same value, 
7 + r (7 includes all relaxation exclusive of bulb loss). Typically, these rates are within a factor of two (see the values 
listed in Fig. |). Then, 
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P'eiax--7P'+ Jl- (14) 

In the bare atom basis, the flux term has a very simple form (Eqn. ^ with no off-diagonal input entries since the 
injected beam has no coherence between the bare atomic states. In the dressed basis, however, there is an injected 
Zeeman coherence, so the flux term takes a considerably more complicated form 

pjiu. = ^F'^ - rp'' (15) 

where F"^ — T^^ (|1)(1| + |2)(2|) T has three diagonal and six off-diagonal entries. 

The Bloch equations are most easily handled by moving to the interaction picture, given by O = g-iHot/h q ^iHot/h^ 
where O is an interaction picture operator. 



C. Steady state solution 



The 4x4 matrix equation (|T^) yields sixteen independent equations that we solve in the steady state. Then, 
the populations in the interaction picture are static, p^^ = 0, and the coherences exhibit sinusoidal precession. In 
particular, 



' , Pib 



R4be^^\ and p4c 



-R4ce**^^+^-'*, where the are time independent, and 
A = LU — UJ24- The other coherences precess at frequencies w^j^ = {Efj_ — E^)/fi. Making these steady state substitutions, 
the sixteen Bloch differential equations transform to a set of time-independent algebraic equations. We assume that 
ujc = W24, SO that the small cavity pulling shift vanishes. The total maser frequency shift is then given by A. 

In terms of dressed basis density matrix elements (rotated out of the interaction picture), the atomic coherence 
Pi2{uj) is given by 



P42(w) 
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(16) 



and the magnetization is found from equation . Inserting this into equation we find the following two conditions 
which determine the maser amplitude and oscillation frequency 



Re [^R,a + SR., - ^i?4c) = -|X24| 



Ith 



(17) 



( -^Ria + SRib - ^^-R4c ) = -\X24\ 



Ith 



where Iq — rVbN is the total atomic flux into the maser bulb and Ith is the threshold flux for maser oscillation with 
our simplified spin-exchange model H: 



Ith 



47r|/i24pQcf7' 



(18) 



Here Vc is the cavity volume and is a dimensionless filling factor 

We numerically solve the time-independent algebraic system of sixteen Bloch equations plus equations ( |l7|) to 
determine the maser frequency shift A as a function of Zeeman detuning 6. We find excellent agreement with the 
previous theoretical bare atom analysis []l2| , within the approximation of equal population decay and decoherence 
rates. (Note that in practice, these decay rates differ by up to a factor of two). 



IV. PHYSICAL INTERPRETATION 



The dressed state analysis provides a straightforward physical interpretation of the double resonance maser fre- 
quency shift. In the absence of the applied Zeeman field, atoms injected in bare state |2) are the sole source of the 
magnetization that provides the positive feedback needed for active oscillation. However, when the near-resonant 
Zeeman field is applied, it also allows atoms injected in the mp = ± 1 states (bare states |1) and |3)) to contribute 
via a two-photon process. A dressed atom interpretation shows how these mp = ± 1 state atoms can become the 
dominant source of maser magnetization as the applied Zeeman field nears resonance. 



5 



Viewed from the dressed atom basis, three factors contribute to this interpretation. First, as shown in Fig. ||(a), 
the apphed Zeeman field shifts the energies of the two dressed levels \a) and |c) symmetrically relative to level \b), 
which remains unperturbed. Second, near the Zeeman resonance, the AF = 1 dipole coupling = (4|/t • Hc|&)^ 
vanishes while and H^^ become equally dominant, as shown in Fig. ||(b). Third, below resonance (S < 0) the 
steady state population of state \a) is greater than that of state |c) {paa > Pec), while above resonance {5 > 0) the 
opposite is true {pcc > Paa), as shown in Fig. ^(c). These dressed state population differences arise from the fact that 
atoms in bare state |1) are injected into the maser while those in bare state |3) are not, under normal operation, so 
in the steady state pn > P33. For large negative Zeeman detunings, \a) — > |1) and |c) |3) (see discussion following 
Eqn. The opposite holds for positive detuning, where \a) — > |3) and |c) 

These three ingredients combine to create the double resonance shift of the maser frequency, shown in Fig. ^(d). 
For small negative Zeeman detunings {\S\ < 2jz), the excess of paa over pcc and the relatively small size of leads 
to maser oscillation primarily on the \a) ^ |4) transition. That is, atoms injected into the maser cavity in the bare 
state |1) contribute significantly to the maser oscillation via a two-photon process: one Zeeman transition photon and 
one microwave photon within the resonant cavity linewidth. This \a) <-> |4) transition is at a slightly higher frequency 
than in the unperturbed (no applied field) maser, so the maser frequency is increased. Conversely, for small positive 
Zeeman detunings (S < 2"fz), the maser oscillates preferentially on the |c) <-*■ |4) transition, and the maser frequency 
is decreased. For larger Zeeman detunings (positive or negative), the coupling of state |4) to unshifted dressed state 
\b) becomes dominant, and the magnitude of the frequency shift is reduced. For zero Zeeman detuning, dressed states 
|a) and |c) are equally populated in the steady state and the maser frequency shift exactly vanishes. 

Injection of an electronic polarization into the maser bulb is needed for the applied Zeeman field to induce a maser 
frequency shift. Since oja and ojc are spaced equally about the unperturbed maser frequency ujt, and since H^^ = 
a necessary condition for a maser shift is a difference in the steady state values of paa and pcc, which is a direct 
consequence of a difference in the injected populations of bare states |1) and |3), i.e., a net electronic polarization. 



V. APPLICATION 



The double resonance hydrogen maser technique was originally studied for use in autotuning the maser cavity p2[ . 
In addition to the double resonance frequency shift, there is a cavity pulling shift for a mistuncd maser cavity, with 
magnitude dependent on the linewidth of the hyperfine transitions, through the line-Q (see Eqn. |^). The applied 
Zeeman radiation depletes the population of bare state |2), thereby increasing the linewidth of the hyperfine transition. 
Andresen showed that the cavity can be tuned to the atomic frequency by modulating the hyperfine linewidth 
with applied Zeeman radiation and adjusting the cavity frequency such that there is no modulation of the maser 
frequency. However, this method requires accurate setting of the applied Zeeman field to the Zeeman resonance (i.e. 
S = 0). 

The double resonance technique can also be used for precision Zeeman spectroscopy in a hydrogen maser. Tra- 
ditionally, the Zeeman frequency in a hydrogen maser operating at low magnetic fields is measured by sweeping an 
audio frequency field through the Zeeman resonance and monitoring the maser power. The power is diminished near 
resonance with a Lorentzian shape with a width on the order of 1 Hz. Typical resolution of the Zeeman frequency 
with this technique is about 100 mHz. However, by utilizing the sharp, antisymmetric profile of the double resonance 
frequency shift, we were able to determine the hydrogen Zeeman frequency with a resolution of about 1 mHz. Recently 
we used this double resonance technique in a search for Lorcntz symmetry violation of the hydrogen atom's electron 
and proton spin motivated by a general extension of the standard model of elementary particle physics [^0|. 



VI. CONCLUSION 



We used the dressed atom formalism to calculate the frequency shift of a hydrogen maser induced by an applied 
field near the F = 1, Amp = ±1 Zeeman transition frequency. The result is in excellent quantitative agreement with 
previous bare atom basis calculations |^,^|, within a simplified spin-exchange approximation and with equal decay 
rates for all populations and coherences. The dressed atom picture provides a simple physical understanding of the 
double resonance frequency shift, and in particular, the atomic polarization dependence of the frequency shift. The 
double resonance technique can be employed in precision spectroscopy of the hydrogen Zeeman frequency, e.g. in a 
test of Lorentz symmetry of the standard model |19,EQ| . 
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FIG. 1. Hydrogen hypcrfinc structure. A hydrogen maser oscillates on the first-order magnetic field-independent |2) ^ |4) 
hyperfine transition near 1420 MHz. The maser typically operates with a static field less than 1 mG. For these low field 
strengths, the two F = 1, Amp = 1 Zeeman frequencies are nearly degenerate, and 1/12 « 1^23 1 kHz. 
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FIG. 2. Hydrogen maser schematic. The solenoid generates a weak static magnetic field Bo which defines a quantization axis 
inside the maser bulb. The microwave cavity field Hp (dashed field lines) and the coherent magnetization M of the atomic 
ensemble form the coupled actively oscillating system. 
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FIG. 3. Double resonance maser frequency shifts. The large open circles are data taken with an input beam of |1) and |2) 
hydrogen atoms. These are compared with Eqn. |^ (full curve) using the parameter values shown. The values of |Xi2| and 
7z were chosen to fit the data, while the remaining parameters were independently measured. The experimental error of each 
measurement is smaller than the circle marking it. The electronic polarization dependence of the double resonance effect is 
illustrated with the dotted data points: with an input beam of |2) and |3) atoms, the shift is inverted. Note that the maser 
frequency shifts for the dotted points were scaled up by a factor of 10 since these data were acquired with a much weaker 
applied Zeeman field. The large variation of the maser frequency shift with Zeeman detuning near resonance, along with the 
excellent maser frequency stability, allows the Zeeman frequency (~ 800 Hz) to be determined to about 3 mHz in a single scan 
of the double resonance such as the dotted data shown here (requiring ~ 20 minutes of data acquisition). 
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FIG. 4. Calculated dressed atom quantities plotted against detuning of the applied Zeeman field (in units of Zeeman linewidth, 
Acjz = 272). The dotted, full, and dashed curves correspond to dressed states |a), \b), and |c), respectively, (a) Dressed atom 
frequencies normalized to the Zeeman Rabi frequency, (b) Interacton Hamiltonian matrix elements (squared) from equation 
(13) in units of (2|/i • Hc|4)^. (c) Steady state populations of dressed states, (d) Fractional double resonance maser frequency 
shift. 
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